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Abstract: We study the phenomenology of a realistic version of the chaotic inflationary 
model, which can be fully and directly explored in particle physics experiments. The 
inflaton mixes with the Standard Model Higgs boson via the scalar potential, and no 
additional scales above the electroweak scale are present in the model. The inflaton-to- 
Higgs coupling is responsible for both reheating in the Early Universe and the inflaton 
production in particle collisions. We find the allowed range of the light inflaton mass, 
270 MeV < m x < 1.8 GeV, and discuss the ways to find the inflaton. The most promising 
are two-body kaon and B-meson decays with branching ratios of orders 10 -9 and 10 -6 , 
respectively. The inflaton is unstable with the lifetime 10~ 9 -10~ 10 s. The inflaton decays 
can be searched for in a beam-target experiment, where, depending on the inflaton mass, 
from several billions to several tenths of millions inflatons can be produced per year with 
modern high-intensity beams. 
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1. Introduction 

In this paper we present an example of how (low energy) particle physics experiments can 
directly probe the inflaton sector (whose dynamics is important at high energies in the 
very Early Universe). 

The common assumption about the inflaton sector is that it is completely decoupled 
from the Standard Model (SM) at energies much lower than the inflationary scale. This 
assumption appears quite natural, since the slow-roll conditions generally permit only tiny 
coupling of the inflaton to any other fields including itself (see e.g. |l|] for a review). This 
implies within the perturbative approach that by integrating out the inflaton sector one 
obtains at low energies some non-renormalisable operators strongly suppressed by both tiny 
couplings and high inflationary scale. These assumptions prevent any direct (laboratory) 
investigation of the inflationary mechanism. 

The situation is quite different if the inflaton sector contains only light fields. In this 
case even weak inflaton coupling to the SM particles can lead to observable signatures in 
laboratory experiments. 

We will concentrate here on the model, which is a particular version of the simple 
chaotic inflation with quartic potential and with inflaton field coupled to the SM Higgs 
boson via a renormalisable operator. We confine ourselves to the case where the Higgs scalar 
potential is scale-free at the tree level, so that its vacuum expectation value is proportional 
to that of the inflaton field and Higgs boson mixes with inflaton. In this model, the inflaton 
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X is found to be light 1 , 270 MeV < m x < 10 GeV, where the lower limit actually comes 
from the searches for decays K — > n + nothing and from the search for axion-like particles 
in CHARM experiment, and the upper limit is related to the requirement of a sufficient 
reheating after inflation. The most promising here are searches of inflaton in decays of K- 
and 5-mesons and searches for inflaton decays in beam-target experiments. 

We present the estimates of meson decay branching ratios into the inflaton, and give the 
estimate for the inflaton production rate in beam-target experiment with beam parameters 
of the T2K, NuMi, CNGS and NuTeV experiments [§, |], gj. The interesting branchings 
start from 10~ 6 for the decay B — > K + x an d 10~ 9 for the decay K — > ir + x ( see 
section ^|). The production rate in the beam-target experiments is of several millions per 
year (at m x ~ 500 MeV) to thousands (at m x ~ 5 GeV). Models with larger masses are 
hard to explore. However, these larger masses correspond to quite heavy SM Higgs boson, 
nih ~ 350-700 GeV. For the described model being consistent up to the inflationary scales 
the Higgs mass should be rrih < 190 GeV, and the inflaton mass is below m x = 1.8 GeV 



(see eq. (2.8)). Hence we conclude that in the model considered here the inflaton sector 
can be fully explored in the particle physics experiments. 

In this respect it is worth mentioning, that hypothetical light bosons, (very) weakly 
coupled to the SM fields, are present in various extensions of the SM (for particular ex- 
amples see ||, ||, [|). Phenomenology of these particles attracted some attention, and in 
Ref. [flCf l the list of relevant experiments is given in the section about searches for axions 
and other very light bosons. These experiments are relevant for the described model (see 
the bounds obtained in this paper). At the same time the theoretical development for our 
case of a light flavour blind scalar dates back to the time, when very light SM Higgs boson 
was still allowed experimentally. Quite a few improvements were done since that time, and 



we mostly used in the current work the results obtained in Refs. [11, 12, 13, 14, [19]. 

The rest of the paper is organised as follows. In section ^ we describe the model and 
outline the viable region in the parameter space, section || is devoted to inflaton decays, 
while meson decays to inflaton are considered in section ||[ There we obtain lower limits 
on the inflaton mass from existing experimental results. In section |B| we study the inflaton 
production in pp-collisions and give production rates for operating high-intensity and high- 
energy beams at JPARC, Fermilab and CERN. Limits on the inflaton mass from the results 
of the CHARM experiment are obtained in section [| There we also give predictions for 
the number of inflatons produced per one year of operation of T2K, NuMi, CNGS and 
NuTeV beams in a beam-target setup. Section [7] contains conclusions. In Appendix [A| 
we discuss consequences of obtained results for a implementation of our model within the 



i/MSM [16], the extension of the SM with three sterile right-handed neutrinos, so that the 
inflaton vacuum expectation value provides sterile neutrino masses. This model is capable 
of explaining neutrino oscillations, dark matter and baryon asymmetry of the Universe 
and, equipped with inflaton sector, provides an example of a full realistic model of particle 
physics. Thus the inflation can be directly tested also in fully realistic extensions of the 
SM. 



1 In fact, there is another window for the inflaton mass, which is also below 1 TeV, see eqs. (2.9), (2.10) 
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2. The model 



We consider the extension of the SM model with an inflaton field introduced in [16]. The 
Lagrangian of this extended model reads 

-CxSM = -CsM + C-XN , 

Cxn = \dpX&>X + l -m\X 2 - f^X 4 - A (rfH - jX 2 )" , (2.1) 

where £sm is the SM Lagrangian without the Higgs potential, while the latter gets modified 
in accordance with (2.1), X is a new neutral scalar field and H is the Higgs doublet. Note, 



that it is supposed that the only scale violating term at tree level is the mass term with 
the negative squared mass — m x for the extra scalar X. This particular choice is sufficient 
to demonstrate the main statement: possibility to test directly the inflationary models in 
particle physics experiments. At the same time, the algebra below is simpler as compared to 
those in the case of general inflaton-Higgs scalar potential, since the number of parameters 
is smaller. 

The SM-like vacuum of the scalar potential in eq. (|2.1j ) is 
(H) = -=, (X) = J-v 



V2 ' N ' V 2a y/20' 
[2a 

with v = J—m x = 246 GeV . (2.2) 
V p A 



The small excitations about this vacuum have the masses 



m h = V2Xv, m x = m h J , (2.3) 



and are rotated with respect to the gauge basis (V2H - v, X) by the angle 

A m v 




(2.4) 



to the leading order in 9 and m x /mh- 

The inflation in the model is supposed to be driven by a flat potential along the 
direction 

H^H^^-X 2 . (2.5) 
A 

The quartic coupling dominates the potential during inflationary and reheating stages. 
Thus, after inflation the Universe expands as at radiation dominated stage, so the number 
of required e-foldings is about N e ~ 62. The normalisation of the matter power spectrum 



1 17] generated during inflation fixes 2 the quartic coupling constant as /3 ~ /3o> where 

A) = 1.5 x 1(T 13 . (2.6) 



2 The estimate ( ]2.6| ) differs from the other values, presented in and in jl^] due to different values of 
iV e used there. 
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In this respect it is worth noting, first, that inflation in this model happens at x ^ Mpi, 
where gravity-induced corrections are expected to be large. We will suppose below, that 
same (yet unknown) mechanism, which guarantees the flatness of scalar potential during 
inflation, operates at x ^ Mpi as well. Thus, the Plank scale physics allows for considering 
the same scalar potential (^4]) at high energies with account of perturbative quantum 
corrections only due to the gauge and Yukawa couplings. These latter are discussed in 
due course. Second, for the quartic inflation there is a tension [17] between predicted 



tensor-to-scalar amplitudes ratio and fits to cosmological data. Though the limits [17] are 



not dramatic yet, weak non-minimal coupling to Ricci scalar, £X 2 R/2, with £ ~ 10 



makes the model fully consistent fl9| , 20, 17] with all current cosmological observations. 



Switching on £ results in a larger value of the quartic coupling constant, /3o < P < 2/?o 
for < £ < 10~ 3 . For our study of inflaton phenomenology at low energies non-minimal 
coupling to gravity is irrelevant, but in further estimates we account for the uncertainty in 
the value of (3, 

P = (1^2) • Pa , (2.7) 

associated with its possible impact. 

Thus, among four parameters in the Lagrangian (|2.1| ), one, P, is fixed by the ampli- 
tudes of primordial perturbations and another combination (|2.2| ) is fixed at the electroweak 
vacuum by the value of the Fermi constant. Two remaining free parameters determine the 
SM Higgs boson mass and the inflaton mass. Further constraints on them are discussed 
below. 

The baryon asymmetry of the Universe is unexplained within the framework of the 
Standard model of particle physics (SM). However, the baryon number is violated at mi- 
croscopic level in primordial plasma, if sphaleron processes are rapid enough (|21| . This 
phenomena is often exploited by mechanisms generating baryon asymmetry within rele- 
vant extensions of the SM. This places a lower bound on the reheating temperature of the 
Universe at the level somewhat above the electroweak scale, and for definiteness we choose 
T r > 150 GeV (for a review see |^]). In the model (|2.l| ), where both quartic coupling 
and Higgs-to-inflaton mixing are very weak, but Higgs boson self-coupling is quite strong 
(A > 0.1 for rrih > 114 GeV), the energy transfer from the inflaton to the SM particles 
is extremely inefficient |l8|]. The stronger is the mixing, the more efficient is the energy 
transfer and the higher is the reheat temperature in the early Universe. Strong mixing 
(larger a) implies lighter inflaton, see eq. ( |2.3| ), and lower bound T T > 150 GeV yields 
upper bounds on the inflaton mass [P 



and 



m h \ / p_ ^ 1/2 

150 GeV/ ' I 1-5 x lO" 13 



m x Z 1-5 • (tTTTT^Tt) • ( i c ) GeV for m x < 2m h , (2.8) 



m < 460 • f — ^ — • ( £ GeV for m y > 2rrih . (2.9) 

x ~ V 150 GeV/ Vl- 5x10 / 

In the latter case the lower bound on the inflaton mass is 

m x > 2m h > 228 GeV . (2.10) 
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In the former case the lower bound on the inflaton mass follows from the upper limit 
on the Higgs- inflaton mixing, appearing from the requirement that quantum corrections, 
originated from this mixing, should not dominate over bare coupling constant /3. With 
the action ( |2.1| ) fixed at the electroweak scale the corrections to the inflationary potential 
/3X 4 /4 can be explicitly calculated and have the form 

m A (X) m 2 (X) 

where m(X) is the mass of the contributing particle in the inflaton background field X 
(taking into account the flat direction (|2.5| ) to obtain the Higgs field background, and \x 
is the electroweak scale. 3 Then, requiring that in the inflationary region X ~ Mp the 
corrections to the quartic coupling j3 are, somewhat arbitrary, smaller than 10%, we get 
from the contribution of the Higgs boson 

a < 1(T 7 , (2.12) 

which precludes large quantum corrections to inflaton quartic coupling driving inflation. 
Limit (|2.12 ) can be converted using ( |2.3[) into the lower bound on the inflaton mass 



i 1 



m h \ 13 V (10 



-7 



2 



^ >u °{t^){tjvw*) {—) MeV - < 2 ' 13) 

Below this limit one should take into account quantum corrections to the inflationary 
potential, which may change the value of the inflaton coupling constant at the electroweak 
scale (^), or even spoil the inflationary picture. 

The proper renormalisation group enhancement of the analysis should be done, once 
any experimental evidence for the light inflaton is found. However, in the larger part of 
the parameter space no significant changes to the described bounds are expected. Limits, 
similar to ( |2.12| ), follow from the requirement of smallness of the SM gauge and Yukawa 
coupling corrections. As far as all the SM particle masses during inflation are proportional 



to y/a/XX (see ( |2.5[ )), all these bounds differ from ( 2.12j ) only by ratios of the coupling 



constants of the form yt/yX, g/V~X, etc. This changes the lower bound ( 2,13|) . Note, that 
the exact value of this bound is not crucial due to the stronger experimental constraints, 
obtained in section |(| However, in some regions of parameter space these corrections may 
lead to much stronger effects — for example, for special ranges of the Higgs boson mass 
the Higgs self-coupling becomes small at inflationary scales, so y/ ' ajXX turns to be large, 
hence so do the SM perturbative corrections to the inflaton potential. 



The weakness of Higgs-to- inflaton mixing ( 2.12j ) is responsible for very tiny inflaton 



interaction with SM particles. This makes searches for heavy inflaton, in the range given 
by ( |2.9D , ( |2.10| ), hopeless in foreseeable future. Indeed, inflaton direct production requires 
to collect and study enormously large statistics in high energy collisions. In contrast, the 
opposite case of light inflaton in the range between (J2j|) and ( gig ) is quite promising, since 



3 Strictly speaking this is the contribution of a scalar particle, while the numeric coefficient changes with 
the number of spin degrees of freedom, and also changes the sign for fermionic particles. 



- 5 - 



inflaton production does not require very high energy at a collision point. Then inflaton 
can be produced in beam-target experiments, where large statistics is achievable. In this 
paper we consider low energy phenomenology of this light inflaton. 

In next Sections we will estimate the decay and production rates of the light inflaton. 
Since inflaton couplings do not depend on the SM Higgs boson mass, its value determines 
only viable inflaton mass range (PI), ( gjg ). Considering for the Higgs boson mass the 



range 114 GeV < rrih < 700 GeV as still possible, in what follows we study the inflaton 
low energy phenomenology for its mass in the interval 

30 MeV < m x < 10 GeV . (2.14) 

Actually, in the model under consideration the upper limit on the SM Higgs boson mass is 
lower than 700 GeV, as we use the same scalar sector (^4|) to describe the inflation at high 
energies. Indeed, if one considers the inflationary model (|2.1|) as it is, it should be valid 



(does not become strongly coupled) up to the energy scale y 1 aj\X ~ 10 15 GeV. As far as 
the inflaton is very weakly coupled with the Higgs field, this requirement is the same, as 
for the Standard Model, leading to the bound < 190 GeV (see, e.g. [pSf] ), or 

m x < 1.8 GeV . (2.15) 



Also, the fit to the electroweak data points at mass interval rrih < 285 GeV ]10| . Neverthe- 
less, in extensions of this model the upper limit on the Higgs boson mass may be higher, 



so we will discuss the whole interval (2.14) to make our study applicable in a more general 
case. 

3. Inflaton decay palette 

Light inflaton decays due to the mixing with the SM Higgs boson ( |2.4| ) . Thus, its branching 



ratios coincide (taking into account the small mixing angle (2.4)) with those of the light 



SM Higgs boson, studied in [11]. In what follows we actually update the results of [11] in 



view of further relevant developments and findings. 

Inflaton of the mass below 900 MeV decays mostly into 

77, e + e", 7r + 7r - , 7r°7r° . 

Mixing (^4|) gives rise to the Yukawa-type inflaton coupling to the SM fermions /, 

C xff = 0^xff = V^ ^Xff , (3.1) 

AJ 1 v m x 

where mj is the fermion mass. Effective inflaton-pion interaction follows from the Higgs 
boson coupling to the trace of the energy-momentum tensor, (cf. fli~4] , [D|) 

C X7T7T = 2k^/2P ■ X- ■ (±d^°& i Tr° + d^ir+d^-) (3.2) 

m X \2 / 

- (3k + 1)V^8 ■ ^ • ml ■ ( ivrV + tt+tt- 
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with m n being the pion mass and 

_ 2N h _ 2 
K ~ 3b "9' 



(3.3) 



where Nh = 3 is the number of heavy flavours, b = 9 is the first coefficient in the QCD 
beta function without heavy quarks. Finally, Higgs-inflaton mixing results in the inflaton 
coupling to VF-boson, which contributes to the triangle one-loop diagram responsible for 
the inflaton decay to a pair of photons. Similar contributions come from fermion loops, so 
that inflaton decay to photons can be described by the effective Lagrangian (cf. p 



^cn^^—xF^, (3.4) 
47r m x 

where a is the fine structure constant, F is a sum of loop contributions from VF-boson and 
fermions / with electric charge eg/, F = Fyy + J2f, colors IfFf 



F w = 2 + 3y[l + (2-y)x 2 ] , (3.5) 
Ff = -2y [1 + (1 - y) : 



ix 2 l 



and y = 4m /mt, with m being the mass of the contributing particle; here also 

x = Arctan , for y > 1 , 

y/y- 1 

x=-\* + i]og 1 _j T — ) ,fory<l. 

In fact, for the interesting range of the inflaton mass the fermion contributions almost 
cancel the VF-boson contribution (the latter dominates over a contribution of each single 
fermion) . 

The inflaton-to-SM fields couplings presented above yield the inflaton decays to leptons 
with the rates: 

r _ = j3 mi mi f 4mj \ 3/2 



x ^ n Air m x \ m| 



(3.6) 



the inflaton decays to pions with the rates: 



_ /3m x (2 llml\ 2 / AmJ 
T x ^- - 2T x ^ nQ + y 1 _ _ , (3.7) 

and the inflaton decay to photons with the rate: 

r^H^) 2 ^. (3.8) 

Note, the tree-level estimate ( |3.7[) is not correct far from the pion threshold, where strong 
final state interactions become important pq| . Thus, in our numerical estimates we follow 
Ref. [15] to improve (T~ 
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If inflaton is heavier than 900 MeV, its hadronic decay modes become more com- 
plicated. First, each quark flavour contributes to the inflaton decay if the inflaton mass 
exceeds the double mass of the lightest hadron of the corresponding flavour. 4 Strange quark 
starts to contribute for the inflaton mass m x > 2rriK, where mx is the iT-meson mass. 
Charm quarks contribute, if m x > 2m_D, where mu is the D-meson mass. Lightest flavour 
states are mesons and close to the thresholds the inflaton decays are described by effective 
interactions similar to (|3.2|). Farther from the thresholds final state interactions become 



important and we follow Ref. [15] to estimate the decay branching rate to pions and kaons. 
This approach is valid until other flavour resonances enter the game or while multimeson 
final states are negligible. For pions this approach becomes unjustified in models with 
inflaton mass above approximately 1 GeV. For mesons of heavier flavours the approach 
fails not far from the corresponding thresholds, where produced mesons are not relativistic 
but many new hadrons enter the game. Far above a quark threshold inflaton decay can be 
described as decay into a pair of quarks due to coupling fl3.1| ) , which subsequently hadro- 
nise. To estimate the decay rate, QCD-corrections should be accounted for, which are the 
same as those for the SM Higgs boson of mass m x . 

Second, if the inflaton mass is close to the mass of some narrow hadronic scalar reso- 
nance (e.g., oniums in cc and bb systems), then interference with such a resonance changes 
significantly inflaton hadronic decay rates in a way quite similar to what was expected for 
the SM Higgs boson if its mass would be close to the bb threshold [24]. We do not consider 



here this case, though the calculations are straightforward. 

Third, for the inflaton mass above 1 GeV the decay into mesons actively proceeds 
via the inflaton-gluon-gluon coupling, which due to the quark one-loop triangle diagram 
similar to that contributing to the inflaton decays to photons. The effective Lagrangian 
describing this decay is [^4|] 



where a s is the strong gauge coupling constant and F = ^2fFf, see (3J5), with sum only 
over quarks. Higher order QCD corrections are also important here, and they coincide with 
those in the case of the SM Higgs boson. The produced gluons hadronise later on. 

Obviously, well above the heavy quark threshold and well above QCD scale the de- 
scription of the inflaton decays in terms of quarks and gluons is well-justified, while in 
the opposite case the effective description in terms of mesons is applicable. Between these 
ranges, where m x ~ 1.5-2.5 GeV, both approximations are not quite correct and no reliable 
description can be presented. At the same time, comparing relevant hadronic contributions 
described within these two approaches valid at lower and upper limits of this "untractable" 
interval, we observe deviations not larger than by an order of magnitude. Thus we con- 
clude, that order-of-magnitude estimates of lifetime and decay rates of leptonic, photonic 
and total hadronic modes for the inflaton mass interval 1.5-2.5 GeV can be obtained by 
some interpolation between these "low-mass" and "high-mass" results. 



4 We do not discuss decays below thresholds with one hadron being off-shell. 
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Hence, for the inflaton decay rates to quarks we obtain the same formula as ( |3.6| ) 
multiplied by the number of colours, 3, and by a factor due to QCD corrections (having in 
mind uncertainties in the value of (3 discussed in section ^, for our estimates we adopt only 
leading order QCD corrections presented in [25|). For decay rates into gluons we obtain 



\F 



2 j Os\ 2 (3 m x 
A-k) Air 



multiplied by the corresponding factor due to leading-order QCD-corrections |25 j. 
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Figure 1: Upper panel: Inflaton branching ratios to various two-body final states as functions of 
inflaton mass m x , for 1.5-2.5 GeV range see discussion in the main text; Lower panel: inflaton 
lifetime t x as a function of the inflaton mass m x . Here we set j3 = fio, while for other values the 
life-time is inversely proportional to f3. With (3 within the favourable interval ( |2.7|) the lifetime can 
be two times smaller. 

The inflaton branching ratios and lifetime are presented in Fig. [l] as functions of the 
inflaton mass. Note, that inflaton partial widths into fermions decrease with increasing 
inflaton mass, while the opposite behaviour is observed for decays into pions due to the 
contribution of kinetic term from the trace of the energy-momentum. Diphoton mode varies 
with the inflaton mass and for the interesting range of parameters reaches its maximum of 
about 2.5% just at muon threshold. Thus, the sub-GeV inflaton predominantly decays into 
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electrons, if its mass is below 200 MeV, otherwise to muons and pions with comparable 
branching ratios. If the innaton mass is below the muon threshold, its life-time is about 
10~ 5 s and above this threshold it falls rapidly down to 10~ 9 s. For the innaton with 
mass well above 1 GeV, heaviest kinematically allowed fermions dominate the innaton 
decay. Decays into photons and electrons are strongly suppressed for m x > 1 GeV, with 
branchings below 10~ 4 . 



4. Inflaton from hadron decays 

In this section we consider the inflaton production in rare decays of mesons and obtain 
corresponding limits on the light inflaton mass. 

First, light inflaton can be produced in two-body meson decays. These are exactly 
the processes widely discussed in the past, when the SM Higgs boson was considered a 



(sub) GeV particle. Making use of the results from [14| one obtains for the amplitudes of 
the kaon decays 



, ,1 3GfV2 

72(1 ~ K) + 



i=c,t ) 



A (K L -> A) = -A (K + -> 7T + X ) , (4.2) 
A {Kg -> xvr ) c e^^^lm £ V^V is , (4.3) 

i=c,t 

with 7i ~ 3.1 x 10~ 7 and much smaller 72. The largest contribution comes from the third 
term in ( [4.11) . At quark level this term is due to the inflaton emission by a virtual quark 
in the quark- W-boson loop. For the branching ratios we get 



Ftotai(^+) lQnM K M k 

,-:-! / 2 IPxl \ Q 2 



1.3 x 10- d • ( ) 9 Z (4.4) 



2 



Br (K L -> n\) » 5.5 x lO" 3 • 



2 \Px\\ ( P\ /100MeV\ 2 



where p x is the inflaton 3-momentum. The branching ratio of K$ is much smaller. For 
the inflaton mass in the kinematically allowed range the squared mixing angle 6 2 is of 



the order 10 5 -10 7 (2.4). For a model with [3 = /3q branching ratios of the kaon decays 



are presented in Fig. together with the relevant existing limits from the searches of the 
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processes K 
inflatons, 



7r + nothing [26 



It follows from Fig. that the models with light 



m x < 120 MeV 



are excluded by negative results of these searches and models with 170 MeV < m x < 
205 MeV are disfavoured. Further increase in sensitivity of these searches by one order 
of magnitude would allow to explore the light inflaton in the mass region 150-250 MeV. 
Thus, kaon decays are the most promising processes to search for the light inflaton. 
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Figure 2: Expected branching ratios Br (K + — > 7r + x), Br {Kl — > n°x) an d experimental bounds 
on Br (K+ -> tt+ X ) from |?J and ]§6). 

In case of a larger inflaton mass heavy mesons have to be considered. The most 
promising here is the rj- meson with the branching ratio of the order (cf. |14], p8[ ) 

while the branching ratios of vector mesons are about two orders of magnitude smaller (cf. 

Two-body decays of charmed mesons similar to K — > ttx are strongly suppressed as 
compared to that decay because of the smallness of the up-quark mass and the correspond- 
ing CKM matrix elements in the amplitude (cf. the third term in ( f4.1|) which dominates). 
Three-body semileptonic decays have larger rates, however they are still quite small |3(| [3l| . 



2 

mr 



«5.7xlO- 9 2 /f^ , 

\ m D J 

f(x) = (1 - 8x + x 2 )(l - x 2 ) - Ylx 2 log(x) . 

Similar formula for the D — > {ivx decay rate is a bit more complicated because of the larger 
muon mass [30], but the rate is suppressed equally strong. On the contrary, decays of the 
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beauty mesons are enhanced as compared to K — > irx- From Ref. [32, ^] and (|2.4j ),(p 
we obtain for the light inflaton: 



2 

2 



\V cb \ 2 \M W ) \ ml) 



(4.7) 



10" 6 



300 MeV\ 



2 




where X s stands for strange meson channel mostly saturated by a sum of pseudoscalar and 
vector kaons. The inflaton with the mass below the muon threshold escapes the detector 
(see Fig. [l]) giving the signatures 

B — ^ K + nothing , B — > K* + nothing . 

Heavier inflaton can decay within the detector, with the most clear mode being the muon 
one at the level of 0.01-1, depending on the mass, see Fig. ||. Note, that collected world 
statistics at b-factories allowed to measure branching fractions of relevant decays B — > 
K^l+l- with accuracy of about 10 7 |&fl| , which is comparable to the expected signal 



Q4.7| ), Thus, an appropriate reanalysis of these data might give a chance to probe the 
inflaton of mass m x ~ 300 MeV. 

Inflaton can also be produced in other meson decays and heavy baryon decays. We 
do not discuss these channels considering them as subdominant for the inflaton production 
and less promising in direct searches for the light inflaton. 



5. Inflaton production in particle collisions 

In this Section we discuss in detail the inflaton production in particle collisions. If the 
collision energy is large enough, the most efficient mechanism of inflaton production is 
kinematically allowed decays of heavy mesons, produced in the collision. Then the produc- 
tion cross section a can be estimated as the product of the meson production cross section 
and the branching of the meson decay into the inflaton, 

-^— = M pp (xs(.0.5BT(K + ^7r + x) + 0.25BT(K L ^TT x)) + XcBr(B^xXs)) , (5.1) 

total 

where total hadron multiplicity M pp and relative parts going into different flavours Xs,c are 
given in Table || for several existing beams, and <7 PPi totai — 40 mbarn is the total proton 
cross section p0|. Here only strange and beauty mesons were taken into account, as they 
give the main contribution. In Fig. |3| we present the estimate of this indirect inflaton 
production in a beam-target experiment for several available proton beams. 

In the models with the inflaton mass above the bottom quark threshold, the dominant 
source of the inflatons is the direct production in hard processes similar to the case of the 
SM Higgs boson. The main channel for the inflaton is the same as for the Higgs boson, i.e. 



the gluon-gluon fusion [38]. The inflaton production is calculated exactly as the production 
of the Higgs boson of the same mass, only it is suppressed by the mixing angle squared 2 . 
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Table 1: Adopted values of relevant for inflaton production parameters of several experiments. 
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Figure 3: The inflaton production cross section in pp collisions normalised to the total pp cross 
section a pp total — 40 mbarn. For all graphs j3 = 1.5 x 10- 13 . The lines are for several existing beam 
energies (see Table 0). For the T2K case the estimate of the direct production by gluon fusion is 
also given, while for higher beam energies it is negligible, compared to the meson decay channel of 
production. At the kaon threshold m x ~ mx the account of the contribution from the 77-meson 



decay (4.6) might somehow smooth the cross section change. 



For several available proton beams the result of the calculation with the parton distribution 
functions from |39| is given in Fig. ||. Here we present the estimate of the direct contribution 
to the inflaton production even for small inflaton mass down to m x > 2 GeV, to illustrate 
the statement that indirect production dominates below the B-meson threshold, cf. Fig. S 

6. Limits from direct searches and predictions for forthcoming experi- 
ments 

The light inflaton will be produced in large amounts in any beam dump, and then decay 
further into a photon or lepton pair, depending on its mass. Search of a penetrating particle 
of this type was performed by the CHARM experiment pCl]. In this article the search for 



a generic axion was performed. However, the only difference between the axion in [40] 
and the light inflaton is in the estimate of the decay rate and production cross section. 
Here we will make rough reanalysis for the case of the light inflaton. As far as in the 
most interesting region of masses the inflaton lifetime is small and a significant amount of 
inflatons decay before they reach the detector, this reanalysis can not be made by simple 
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Figure 4: The ratio of the inflaton direct production cross section in pp collisions to the total pp 
cross section, a pp to tai ~ 40 mbarn, for inflationary j3 — /3q. 



rescaling of the result. So, first we will reproduce in the simplest possible way the resulting 
bound from [f5|, and then change in the analysis the decay and production rates. The 



resulting picture in [40[ can be reproduced by demanding that the number of decays in the 
detector is larger than the background. The number of decays in the detector is roughly 
estimated as 



N ~ N Q a x e~ 



1 



(6.1) 



where Nq is the overall coefficient describing luminosity, ax is the production cross section 



of the axion (eq. (5) in [40]), T is the decay width (sum of eqs. (3) and (4) for muons and 
electrons in |40]]), Idee = 480 m is the decay length before the detector, Idetector = 35 m is 
the detector length, 7 = E/mx is the typical relativistic gamma factor of the axion, with 

10- 17 . 



E 



■> 10 GeV. Then, Fig. 4 from [40] is approximately reproduced for N / [Nqg^o) 
Using the same logic we can obtain the bound for the inflaton. We then get instead of 



N ~ N a prod (Br{ X -> 77) + Br(x -> ee) + Br(x -)■ up)) 



p L dec 

x e c t 



1 



, (6.2) 



where for the inflaton production cross section a x we take (5.1), and we adopt the simple 
estimate of the ttq yield: a n o /o" pPj totai = M pp /3, with M pp = 13 for the CNGS beam (see 
Table [l]). Then, the region forbidden by the experiment |4C] is given in Fig. |5|. 

One can see, that for the inflationary self coupling (3 = /3o the CHARM results exclude 
the inflaton with masses m x > 280 MeV, while for the upper limit of the interval ( |2,7| ) we 
have 

m x > 270 MeV . (6.3) 

One should be somewhat careful taking this limit. The exact dependence of the bound 
on /3 in this region is sensitive to the proper analysis of the CHARM experiment (proper 
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Figure 5: Limits on the inflaton properties from the CHARM axion searches. The region of (m x , (3) 
surrounded by the curve is forbidden. One can see, that for inflationary f3 — /3 the inflatons with 
masses up to approximately 280 MeV are ruled out. 

simulation of energy distribution of the produced inflatons, detector sensitivity to different 
decay modes, etc.). From the plot in Fig. [| we conclude, however, that the value of 
m x > 210 MeV is a conservative bound: within the simple approach above we get an order 
of magnitude "safety margin". The definite conclusion about higher values of masses for 
/3 ~ /3o requires careful reanalysis of the CHARM data. 

For the allowed inflaton mass range and using the results obtained in section [| we esti- 
mate the number of inflatons produced during one year of running at designed luminosity 
for several experiments by multiplying the cross section ratio by the number of protons on 
target, see Fig. |6| and Table |l[ This number should be taken with a grain of salt, as far as 
we have totally neglected all possible geometrical factors and possible collimation or deflec- 
tion of the produced charged kaons. However one can conclude that higher energy beams 
are certainly preferable in the searches for the inflaton, because the dominant production 
mode is decays of beauty hadrons. 5 For small mass the number of the inflatons per year 
can exceed several millions, while at m x ~ 5 GeV it is about a thousand at best. 

7. Conclusions 

In this paper we presented an example of a simple inflationary model, which can be fully 
explored in particle physics experiments. The inflaton is light and can be hunted for in the 
beam-target experiments and in two-body meson decays. The inflaton life-time is about 
10~ 9 -10~ 10 s, so it decays close to the production point. The dominant mechanism of 
inflaton production is the meson decays. Thus, and it is quite amusing in fact, the first 
stage of the Universe evolution described by means of quantum field theory can be directly 
tested with modern experimental techniques. 

5 This is not so for low mass inflatons, which can be produced in kaon decays and one can win slightly 
because of usually higher luminosity of the lower energy beam. 
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Figure 6: Number of the inflatons produced with one-year statistics for the reference beams. All 
geometric factors are neglected. The vertical line indicates the lower limit on the inflaton mass 
protect (6.3). 



The mass limits on the light inflaton are found to be 270 MeV < m x < 1.8 GeV. Note, 
that careful reanalysis of the results of the CHARM experiment [^(J may change the lower 
bound given here. 

The paper analysed the quartic inflaton self-interaction with minimal (or very weakly 
non-minimal) coupling to gravity, which is still allowed by the current WMAP data. In near 
future results from the PLANCK experiment will determine the inflationary parameters — 
spectral index n s and tensor-to-scalar ratio r — with much better precision. This will allow 
to fix the value of the non-minimal coupling constant £ (see relations between £ and n s , r 
in |l8|, |(], |j~9[|). In case it is large, the laboratory detection will still be possible, while the 
analysis of the current work should be rescaled for larger value of the self-coupling constant 
f3 and thus for stronger signal. 



Note also that the inflaton model discussed in this paper has been suggested in Ref. [16] 



as an extension of the z^MSM [41, |42j. In that joint model the inflaton vacuum expectation 
value provides the three sterile neutrino of z/MSM with Majorana masses, while decays of 
inflatons in the early Universe to the lightest sterile neutrino produce the dark matter. In 
Appendix [A| we check, that the limits presented in z/MSM itself do not limit the interval for 
the inflaton mass obtained in this paper. On the contrary, the upper limit on the inflaton 
mass constrains the mass of dark matter sterile neutrino in this model. This proves that 
inflation can be directly tested in fully realistic extensions of the SM. 

We are indebted to M. Shaposhnikov for valuable and inspirational discussions through- 
out the work, P. Pakhlov for discussion on recent Belle results, V. Rubakov for valuable 
comments. D.G. thanks EPFL for hospitality. The work of D.G. was supported in part 
by the Russian Foundation for Basic Research (grants 08-02-00473a and 07-02-00820a), 
by the grants of the President of the Russian Federation MK-1957.2008.2, NS-1616.2008.2 
(government contract 02.740.11.0244), by FAE program (government contract II520), by 
the Scopes grant of Swiss National Science Foundation and by the grant of the Russian 
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A. The z/MSM extension 



In Ref. [16] the inflaton model considered in this paper has been implemented in the 



framework of vMSM [41, 42 1. 



Neutrino oscillations is the only direct evidence of physics beyond the Standard Model 
of particle physics. Cosmology provides two other evidences — dark matter and baryon 
asymmetry of the Universe — which are not explained within SM if General Relativity is 
a correct theory of gravity. These three problems can be addressed within a neutrino 
minimal Standard Model (i/IVISM) J4l| , B^ ], suggested as a minimal version of SM capable 
of explaining all of the three. Enlarged by the additional scalar field |l(| coupled to the 
SM Higgs boson via tree-level-scale-invariant interaction, z^MSM provides both early-time 
inflation and common source of electroweak symmetry breaking and of sterile neutrino 
masses. This one-energy scale model 6 is a minimal, full and self-consistent extension of the 
Standard Model of particle physics. 

This extension implies addition to the model Q2.1| ) three new singlet right-handed 
neutrinos Ni, I = 1, 2, 3 with the Lagrangian 

£,msm = iNil^Nj + ^F aI L a NjH - ^NfNjX + h.c. 

where L a (a = 1, 2, 3) are charged lepton doublets and H = eH* , where e is 2 x 2 antisym- 
metric matrix and H is the Higgs doublet. 

To begin with, let us discuss the allowed ranges of the Yukawa coupling constants // 
(Yukawas Fj a are irrelevant for this study). The flatness of the inflaton potential implies 
smallness of the quantum corrections to the quartic coupling Q2,6| ) . Requiring again smaller 
than 10% contribution we obtain 

// < 1.5 x 1(T 3 (A.l) 
and get for the sterile neutrino masses 

1.5xl0- 13 \ 1/2 / fj 



Mi = fi (X) < 270 ■ * — - ? GeV . (A.2) 

J x 1 VlOOMeV/V P J V 1 - 5 * 10 / 

A lower limit on the masses of the two heavier sterile neutrinos, -^2,3, follows from suc- 
cessfulness of the Big Bang Nucleosynthesis (BBN), that does not constrain the inflaton 
mass. 

The lightest sterile neutrino in the model ( |2.1[ ) can be stable at cosmological time-scale 
and comprise the dark matter of the Universe. In this case, the sterile neutrino is quite 
light, Mi < 1 GeV, and should not get thermalised in primordial plasma, to be viable dark 
matter. The latter is natural, since stability at cosmological time-scale implies extremely 
small (if any) mixing with active neutrino via F\ a , and small mass implies tiny coupling 



'For a discussion of gauge hierarchy problem in models of this type see, e.g. Ref. 
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to inflaton. Light inflatons are in equilibrium in the primordial plasma and decay to the 
lightest sterile neutrino mostly at the temperature T ~ m x . They provide contribution of 
the sterile neutrino N% to the energy density of the Universe today 

1.6/(m x ) ( Mi \ 3 /100MeV\ 3 

^ = " 1.5 x 10-13 ) '\~^~) ' ^ 

where S > 1 is a dilution factor accounting for a possible entropy production due to the 
late decays of the heavier sterile neutrinos [44| and the function / (m x ) is determined by 



the number of degrees of freedom (T) in a primordial plasma at inflaton decays. It 
changes monotonically from 0.9 to 0.4 for inflaton mass from 70 MeV to 500 MeV and for 
heavier inflaton can be approximated as / (m x ) ~ [10.75/g* (m x /3)] 3 ^ 2 . Other mechanisms 
can also contribute to the dark matter production. Hence, at a given m x equation ( |A.3| ) 
implies the upper limit 

M ' <(^Y 3 (m 1/3 ( i -^-) 1/3 (A)- (A.,, 



10keV \6Af(m x )J \0.25 J \ p J VlOOMeV 

The corresponding upper limit on /i supersedes (|A.1|) . 

There are several notes in order. First, from ( |A.4| ) one concludes, that with the inflaton 
mass in the range (|2.§|), ( [2.13 ), the lightest sterile neutrino can naturally be a Warm Dark 



Matter candidate. Second, with model parameters tuned within their ranges to minimise 



r.h.s. of the inequality (A.4), one does not exceed the lower bound on the sterile neutrino 



mass, Mi > 1.7 keV, from the study [45, [5| of the dark matter phase space density in 
dwarf spheroidal galaxies. Hence, in this model limits from the dark matter sector do not 
shrink the allowed parameter region, and in particular, the inflaton mass range. Third, 
maximising r.h.s. of inequality (A.4) one obtains an upper limit on the lightest sterile 
neutrino mass, 

Hence, given the range Q2.15| ) the lightest sterile neutrino is significantly lighter than elec- 
tron (and other charged SM fermions). Since the light inflaton decays to SM fermions due 
to mixing with the Higgs boson, these partial decay rates are proportional to the squared 
masses of the corresponding fermions similar to its decay rates to sterile neutrinos. The 
lightness of the dark matter sterile neutrino guarantees, that inflaton decays to the dark 
matter neutrino never suppresses its decay branching ratios to visible channels. 

From the formulas above one concludes that the limits on the inflaton mass obtained 
in this paper are not affected due to additional constraints typical for ^MSM. At the same 
time, the limits ( |A.5[ ), ( p. 15 ) imply that the dark matter is lighter than about 100 keV and 



with account of the allowed Higgs boson mass (see Sec. ||). Once the inflaton is found the 



upper limit ( A. 5 ) on the dark matter neutrino mass will be settled. And vice versa, once 
the dark matter neutrinos are found, eq. ( |A.5| ) fixes a lower limit on the inflaton mass; this 
limit can supersede limit (|6.3| ) from direct searches, if Mi > 10 keV. 
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Note that inflaton can decay into sterile neutrinos ^2,3 if kinematically allowed. The 
decay mode into the not-dark-matter sterile neutrinos (invisible mode) has the width 

This formula is very similar to ( |3.6|) . This decay mode can even be dominant, if the sterile 
neutrinos A^ are the heaviest fermions between the kinematically allowed ones. As far 
as M2,3 < rn-n is disfavoured by the BBN, this can be relevant for inflaton masses above 
approximately 300 MeV. As a consequence, inflaton lifetime can be somewhat shortened 
as compared to the results presented in Fig. |l]b. 
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